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I.  INTRODUCTION 


Beginning  October  1,  1978  a  research  grant  (AF-AFOSR-79-0009)  was 
awarded  to  the  Polytechnic  Institute  of  New  York  by  AFOSR  with  Professor 
B.R.  Cheo  as  the  Principal  Investigator.  Since  October  1,  1979,  Profes¬ 
sor  S.P.  Kuo  shared  responsibilities  as  the  Co-Principal  Investigator  and 
became  Principal  Investigator  in  October  1982.  Professor  B.R.  Cheo  then 
assumed  responsibilities  as  Director  of  the  Laboratory  of  Plasma  Engineer¬ 
ing  and  Science  as  well  as  Co-Principal  Investigator  on  the  subsequent 

separate  OSR  program  in  mm  wave  generation  beginning  October  1,  1982. 

w  £ 

Essentially  five  topics  been  studied  during  this  time  period;  three 

are  essentially  experimental  investigations,  and  the  other  two  are  theo¬ 
retical  in  nature.  The  five  tonics  include: 

1)  Propagation  of  microwaves  along  a  plasma  column  and  harmonic 
generation  of  electrostatic  ion  cyclotron  waves, 

2)  Imploding  tube  experiment, 

3)  RF  Generated  Current  in  a  Magnetized  Plasma  Using  a  Slow  Wave 
Structure, 

4)  Wave-particle  interaction  at  cyclotron  resonances,  and 

5)  Turbulent  interaction  between  waves  and  charged  particles.  , 

Topics  1-3  comprise  the  experimental  investigation  and  the  last  two 
are  theoretical.  In  the  next  section  we  summarize  the  achievements  made 
on  each  of  these  topics. 


II.  PROGRESS  REPORT 

1)  Harmonic  Generation:  Two  independent  investigations  have  been  con¬ 
ducted  in  Polytechnic  Hollow  Cathode  Discharge  plasma  device.  In 
the  first  investigation  the  phenomena  of  harmonic  generation  of  the 
electrostatic  ion-cyclotron  wave  and  suppression  of  the  second  and 
enhancement  of  the  third  harmonic  wave  have  been  observed  and  a 
theoretical  model  to  explain  it  is  proposed.  The  agreement  between 
the  two  is  favorable. 

In  the  experiment  the  electrostatic  ion-cyclotron  wave  is  ex¬ 
cited  parametrically  with  the  second  harmonic  electron-cyclotron 
wave  by  injecting  a  monochromatic  electromagnetic  wave  at  fQ=9.23  Ghz. 
The  generation  of  the  second  and  of  the  third  harmonic  waves  is  clear¬ 
ly  observed  on  the  spectrum  analyser.  The  evolution  process  of  sup¬ 
pression  of  the  second  harmonic  and  enhancement  of  the  third  har¬ 
monic  wave  is  studied  by  varying  a  d.c.  voltage  from  zero  value  on 
the  probe  that  is  used  to  pick  up  the  signal. 

Starting  with  the  Maxwell's  equations  for  electrostatic  wave  we 
obtained  the  rate  equations  for  the  second  and  the  third  harmonics. 

The  nonlinear  coupling  coefficients  can  then  be  calculated  by  using 
the  Vlasov  equation.  After  a  lengthy  manipulation  we  examine  the 
steady  state  solutions  of  the  rate  equations  which  are  compared  favor¬ 
ably  with  the  experiments. 

In  the  second  investigation,  by  adjusting  the  operating  para¬ 
meters  of  the  device,  two  distinct  modes  of  plasma  column  can  be  ob¬ 
tained.  The  first  one  is  high  density  (a>  <  u)  )  with  bigger  radius 

pe  o 

and  a  brighter  periphery  than  the  interior.  The  second  one  is  of  low 
density  (w  <  wq)  with  a  small  radius  and  the  apparent  light  intensity 
is  maximum  at  the  center  of  the  column.  The  microwave  intensity  is 
maximum  on  the  surface  of  the  plasma  column  for  the  first  mode.  The 
microwave  intensity  is  maximum  at  the  center  of  the  column  for  the 
second  mode. 
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A  simple  theoretical  model  has  been  proposed.  The  plasma  is 
considered  as  cylindrical  dielectric  wave  medium  which  is  situated 
in  a  longitudinally  applied  d.c.  magnetic  field.  Full  Maxwell's 
equations  have  been  used  to  examine  the  problem  instead  of  the  usual 
simplified  electrostatic  approximation.  Agreement  between  theory 
and  observations  is  favorable  for  both  modes. 

Project  is  completed.  A  technical  report  is  issued:  (H.M. 

Huang  and  B.R.  Cheo,  Propagation  and  Nonlinear  Studies  of  Electro¬ 
magnetic  Wave  in  a  Self-Sustained  Plamsa  Column,”  Scientific  Re¬ 
port  No.  POLY-EE-80-005 ,  June,  1980). 

2)  Imploding  Tube  Experiment:  An  investigation  of  electrically  explod¬ 
ing  metallic  thin  shells  is  conducted  experimentally  and  theoretically. 
An  instability  is  identified  experimentally  on  the  tube  prior  to  vapor¬ 
ization.  A  theoretical  analysis  is  developed  that  predicts  this  in¬ 
stability. 

The  metallic  shells  are  exploded  via  a  fast  capacitor  bank.  The 
tube  is  observed  with  a  modified  shadowgraph  technique  and  a  long  wave¬ 
length  (=2.3  cm)  instability  is  observed  to  dominate  at  the  onset. 
Theoretical  treatment  of  the  problem  is  derived  in  part  from  the  Z- 
pinch  instability  method  of  Kruskal  and  Tuck  and  the  mechanical  stress 
analysis  is  adapted  mainly  from  the  techniques  of  Onat  on  thin  shelled 
tubes.  These  methods  are  integrated  to  form  a  theoretical  basis  for 
the  instability  analysis.  The  experimental  observations  and  the  re¬ 
sults  of  the  theoretical  analysis  compare  favorably. 

Project  is  completed.  A  technical  report  is  issued:  (Charles 
D.  Hechtman  and  B.R.  Cheo,  "Investigation  of  Electrically  Exploding 
Thin  Shelled  Metallic  Tubes,"  Scientific  Report  No.  POLY-EE-80-007, 
October,  1980). 


3)  RF  Driven  Currents:  The  work  in  this  section  has  not  been  reported 
previously.  This  experimental  investigation  has  reached  conclu¬ 
sion  and  a  scientific  report  (POLY-EE-83-004)  is  being  prepared  and 
will  be  issued  shortly.  This  work  was  motivated  by  the  possible 
steady-state  operation  of  plasma  devices  by  using  unidirectional  rf 
waves  to  drive  dc  currents.  Experiments  were  performed  in  our 
linear  ECRH  plasma  device  which  uses  a  set  of  phased  capacitor  plates 
excited  by  a  3.25  MHz  rf  source  to  excite  unidirectional  electro¬ 
static  (w/k  <<  c,  ft.  <<  w  <<  ,  ft  )  waves.  Fig.  1  shows  the  basic 

z  i  pe  e 

experimental  setup  used.  The  slow  wave  of  structure  is  located  in 
the  region  of  uniform  magnetic  field  and  is  located  at  the  low  dens¬ 
ity  periphery  of  the  plasma  column.  Fig. 2  shows  additional  detail 
of  the  experimental  arrangement.  Diagnostics  used  in  the  experiment 
work  include: 

1)  Longitudinally  and  radially  movable  Langmuir  probes  used  to  deter¬ 
mine  electron  temperature,  electron  density,  and  space  potential. 

2)  Radially  and  azlmuthally  movable  current  probes,  used  to  measure 
radial  and  azimuthal  profiles  of  current  density  in  the  z-direc- 
tion. 

3)  Electrostatic  Energy  Analyzer,  used  to  measure  the  parallel 
electron  energy  distribution  function. 

4)  RF  potential  probe,  used  in  conjunction  with  rf  interfero¬ 
meter  to  obtain  amplitude  and  phase  characteristics  of  the  rf 
DOtential  associated  with  electrostatic  waves  in  the  plasma. 

It  is  found  experimentally  that  the  electron  flow  is  localized  rad¬ 
ially  and  azlmuthally  to  the  lower  hybrid  resonance  layer  near  the  outer 
region  of  the  plasma  column.  Large  modifications  in  the  electron  energy 
distribution  function  are  also  found  in  this  region  as  well  as  large  radial 
phase  variations  in  the  rf  potential  in  this  region.  These  modifications 
depend  strongly  on  the  sign  of  k^.  Detailed  experimental  results  and  their 
interpretation  will  be  presented  in  a  scientific  report  (POLY-EE-83-004) 
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(Ph.D.  Dissertation,  Brian  R.  Poole)  to  be  issued  in  the  early  Fall  of 
1983.  Results  of  the  experimental  investigation  were  presented  at  the 
1983  IEEE  International  Conference  on  Plasma  Science. 
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FIG.2  DETAILED  EXPERIMENTAL  SETUP 


CURRENT  SENSING  RESISTORS 


4)  Wave-Particle  Interaction  at  Cyclotron  Resonances: 

a)  Perpendicular  Electron  Cyclotron  Resonance  Heating  by  Ordinary 
Modes (2) 

Electron  heating  at  resonance  by  the  ordinary  mode  is  analyzed. 
The  equations  of  motion  are  integrated  along  the  resonance  tra¬ 
jectory  showing  no  parallel  heating.  Evolution  of  the  distribu- 

2 

tion  function  shows  that  grows  expoentially  in  t  for  second 
harmonic  and  algebraically  at  the  fundamental. 

b)  The  Effect  of  the  Bouncing  Motion  of  Electrons  on  Electron 
Cyclotron  Resonance  Heating^ 

Cyclotron  resonance  heating  of  bouncing  electrons  by  an  ob¬ 
liquely  incident  wave  field  is  analyzed.  Continuous  interac¬ 
tion  between  electrons  and  wave  fields  throughout  the  entire 
range  of  electron  excursion  has  been  considered  in  the  analysis. 
The  results  show  that  the  bouncing  motion  of  electrons  serves  to 
alleviate  the  detuning  effect  of  frequency  mismatch  and  efficient 

heating  can  be  achieved.  Kinetic  temperature  is  found  to  in- 

2 

crease  algebraically  in  t  for  the  fundamental  resonance  and 
ej  aentially  in  t  for  the  second  harmonic  heating.  Due  to  the 
interaction  between  electrons  and  the  parallel  component  of  the 
wave  electric  field,  heated  electrons  tend  to  focus  to  the  mid¬ 
plane  and  their  excursion  amplitudes  also  become  filamented  in 
the  steady  state. 

Copies  of  these  journal  articles  can  be  found  in  the  appendix 
of  this  report. 

5)  Turbulent  Interaction  Between  Waves  and  Charged  Particles: 

Turbulent  Heating  of  Parametric  Instabilities  in  Unmagnetized 
(3) 

Plasmas 

Consideration  of  the  effect  of  a  uniform  pump  field  on  the  particle 
orbits  in  a  Vlasov  plasma  leads  to  a  modified  diffusion  coefficient. 
When  the  particles  oscillate  in  the  pump  field,  the  turbulent  wave 
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phase  velocity  seen  by  the  particles  is  Doppler  shifted  by 
multiples  of  u>o/k  is  the  pump  frequency).  Hence,  strong 
interactions  between  the  particles  and  various  components  of 

the  turbulent  field  will  take  place.  It  is  shown  that  when 

1/2 

the  pump  field  is  sufficiently  strong  Eq  >  (Ann^^T^)  ,  bulk 

heating  can  dominate  over  tail  heating  and  excitation  of  elec¬ 
trons  to  energy  levels  higher  than  the  normal  suprathermal  values 
is  possible.  This  field  strength  is  within  the  range  of  laser 
fusion 
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This  section  contains  important  journal  publications  relating  to 
the  work  on  cyclotron  resonance  and  turbulence. 
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Perpendicular  electron  cyclotron  resonance  heating  by 
ordinary  modes 

S.  P.  Kuo  and  B.  R.  Cheo 

Polytechnic  Institute  of  Sew  York.  Formintdale.  Sew  York  1 1 735 
(Received  6  May  1980;  accepted  8  January  1981) 

Electron  heating  at  resonance  by  the  ordinary  mode  is  analyzed.  The  equations  of  motion  arc  integrated  along 
the  resonance  trajectory  showing  no  parallel  heating.  Evolution  of  the  distribution  functions  shows  that  Tt 
grows  exponentially  in  t’  for  second  harmonic  heating  and  algebraically  at  the  fundamental. 


Recently,  there  has  been  renewed  interest  in  electron 
cyclotron  resonance  heating1  **  of  plasmas  in  connection 
with  the  need  (or  electron  heating  in  magnetically  con- 
fined  devices,  and  with  ionospheric  heating  experi¬ 
ments.  It  plays  a  crucial  role  in  the  operation  of  the 
Elmo  Bumpy  Torus. It  was  thought  that  because  of 
the  parallel  polarization  of  the  electric  field  of  the  heat¬ 
ing  wave,  ordinary  modes  would  only  lead  to  parallel 
heating.  However,  in  the  following,  we  shall  show  that 
at  resonance,  when  the  effect  of  the  oscillating  mag¬ 
netic  field  is  included,  there  is  strong  perpendicular 
heating  and  no  parallel  heating  at  all. 


v'U)  =  u,o  exp[-j(fW  +  0O)] — — Lexp(-iflf) 

x  j  ds  u,(s)exp(tns)cos[fcr(s)-  us],  (4) 

v-{t)  =  v«+9t  l'ds(v-{s)-Tt) 

Xcos[fex(s)~  ws] ,  (5) 

where  f X  =  eBt/mc  and  the  initial  conditions  r,  =  r(0), 
v,  =  v(0)  =  ul0(x  cosS0  +  y  sine,)  +  (0) , 

and 


We  follow  standard  procedures  to  determine  particle 
motion  in  the  heating  wave  field  of  an  ordinary  mode  in¬ 
cluding  the  effect  of  the  oscillating  magnetic  field  of  the 
heating  wave  in  the  Lorentz  force.  A  perturbed  (self- 
consistent)  trajectory  is  found  and  used  in  the  integra¬ 
tion  of  the  equations  of  motion  at  resonance  (<j  =nfi). 

It  is  found  that  the  parallel  heating  term  vanishes  ex¬ 
actly  at  resonance  in  agreement  with  previous  work1*’11 
using  a  first-order  quasi-linear  diffusion  theory  (weak 
heating),  but  differs  from  the  result  of  single  particle 
theory*  which  neglected  the  oscillating  magnetic  field. 
Efficient  perpendicular  heating  is  predicted.  The  tem¬ 
poral  evolution  of  the  electron  distribution  function  and 
of  the  perpendicular  temperature  are  derived  for  the 
fundamental  resonance  and  for  the  second  harmonic. 

In  a  homogeneous  plasma  background  imbedded  in  a 
uniform  field  Bs,  and  an  electromagnetic  field  having 
ordinary  polarization,  the  electron  motion  obeys: 


dr{t) 

dt 


v(t) , 


(1) 


H  (t)  =  (ix  +  yy)  sinil/  -  try  -  yx)(l  -  cos  fit)  +  ziflt . 

In  order  to  find  the  exact  trajectory,  (3)— (5)  must  be 
solved  simultaneously.  However,  at  resonances  u>  =  nft, 
the  response  to  the  resonance  interaction  dominates  all 
other  terms.  The  exact  trajectory  would  then  be  that 
given  by  the  resonant  (secular)  term  modulated  by  ne¬ 
gligibly  small  oscillations  due  to  the  other  terms  which 
are  smaller  by  a  factor  (fl/)"*.  In  this  limit,  we  have: 

x(/)=x,  +  [if (/)/*] sin9(/)-  (^,/n)  sinS, ,  (6) 

where 

w(f)  =  [(i?  +  M)’  +  (fea),],/’,  »(/)  =  9,  +  nz  +  *>(/),  (7) 

«(/)  f  d/'{/,.i(w')sin[*  +  (n-lWi'l 

z»0  ■'0 

-  J.,,(if')sin[*  +  (n  +  l)i(i']},  (8) 

8(/)  =  -^dL  dt'{j,.\{u>')  cos[*  +  (n  -  l)d>'} 

Z&o  •'0 

+  ./„,,(«(')  cos[*  +  (n +  !)<(>']},  (9) 


=  -^(E[r(/),/]+^x{B[r(/),/]  +  B,}),  (2) 

where  E(r,  ()  and  B(r,f)  are  oscillating  fields  of  the 
wave.  We  choose  a  coordinate  system  in  which 

B,  =  £fi, ,  k  =  xk,  E(r,  /)  =  *E|  cos(kx  -  w/) , 

and  B(r,  /)  =  -yB,  cos(fcx  -  w /),  where  Bt  =  kcEx/u.  Let 
v~=v,~  ivf,  hence,  (1)  and  (2)  can  be  integrated  as 

r(/)  =  r,  +  H^‘v»  — ^  j £  ds  «(/-*)•  ^E(r(s),  s] 

+  ^XB[r(s),*|^,  (3) 


with  *  =*x,-Tjsin 80  +  n60,  i)*kv0/n,  iv'=w{t'),  = 

and  4>(t)  =  tan*'[ka/(i)  +  *£))].  Substituting  (6)  and  the 
resonance  response  part  of  v,  =  [n«//*)cos8(/)J  into  (5), 
yields 

t '.(f)  =  v.,  +  ^exp(i> )  dt'(vl%  +  n  0(f) 

x  {J,.|(w')exp[«0«  +  l)tf>'  J  + 
x  exp[t(n  -  1)*']}-  -1  «/„<w')exp(«»*') 

°i 

-  ifla (/')(«/..| (if ')exp[i(n  +  l)*']  -  J,.x{w') 

xexp[t(n-  l)f']^exp[-i(w-»in)/']+c.c.^.  (10) 
If  the  oscillating  magnetic  field  is  neglected  and  the  un- 
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perturbed  trajectory  is  used,  previous  results  can  be 
recovered.*  Using  the  trajectory  (6),  Eq.  (10)  becomes 

v,(/)  =  v,%  +  »»  +  *W0')l 

*{  cos[*  +  (n  + 

+  coa[*  +  0t-l)r]J.H(*')}- (2»n/k)  cos(e  +  «#') 

X  +  fta(/'){sln[*  +  (n  +  l)^']J.M(w') 

-  sin(«  +0t- 

=  £  df  w‘  cos(*  +  it#')[</,«i(u>') 

+  </„.((*>')  -  (2n/ui'y.(u>')]  =  vH .  (11) 

Here, the  relations  cos# =  (q  ♦  t0)/w ,  sin#  *  ka/w,u>  «nO, 
and  the  Identity  Jm^(u>)  +  -  (2n/w)J,(w)  =  0  have 

been  employed.  It  is  seen  that  there  is  no  parallel 
resonant  heating  at  all.  Physically,  this  may  be  ex¬ 
plained  as  follows:  When  an  electron  is  moving  in  reso¬ 
nance  with  the  field  along  the  resonant  trajectory,  it 
sees  an  electric  field  consisting  of  a  constant  term  plus 
a  series  of  oscillations  at  nfl.  The  constant  term  of 
the  field  is  proportional  to  «/„( u>).  In  addition,  in  the  s 
direction  the  electron  also  experiences  a  vx  B  force 
consisting  of  two  constant  terms  plus  oscillations.  The 
two  constant  terms,  one  upshifted  and  one  downshifted, 
proportional  to  respectively,  would  exactly 

cancel  that  due  to  the  electric  field.  The  net  Lorentz 
force  on  the  electron  in  the  parallel  direction  is  hence 
oscillatory  and  does  not  contribute  to  heating. 

Substituting  (11)  into  (4),  we  further  obtain 
o‘(f)  =  [nw(t)/fc]exp[-»0(t)],  vj(f)s*(nx>/k)*.  (12) 

With  the  aid  of  (8)  and  (9)  and  the  trajectories  (6),  (1), 
(11),  and  (12),  the  temporal  evolution  of  the  perpendicu¬ 
lar  energy  of  electrons  and  the  electron  distribution 
function  can  be  determined.  In  the  collisionless  limit, 
the  evolution  of  the  distribution  function  follows  the 
equation  d/,(v,r,t)/dt  =  0  along  the  trajectories  (6),  (7), 
(11),  and  (12).  This  equation  gives  the  relation 
f.(v,  r,  f)  =/,(▼»,  r»,  0)  provided  that  (v,  r)  are  related  to 
(v»,  r»)  through  the  relations  (6),  (7),  (11),  and  (12). 
Thus,  for  a  given  initial  distribution  f(vt,  r«),  we  may 
determine  the  time  evolved  distribution  function  under 
the  resonance  heating  limit.  The  kinetic  perpendicular 
temperature  Tx{t)  is  thus: 

TAt)  ■  J  f  +  n0O)F + (nor  (r)I*}f.(yo.  r«,  o)</vc .  (13) 

We  consider  two  cases  for  an  initial  Maxwellian 
/,(v„  ro,0)  -  frn/2»T0)v*  exp(- W27'0)(u*X0+  01: 

(i)  Fundamental  cyclotron  resonance  heating  (n  =  l): 
Let  «/,(*>)*  1  and  J,(ie)*  0  for  lwl«l,  (wav,/c),  thus 
aU)  a  lB,v*i/2B,)  sins  and  0(f)  =  - (B,urtf/2Be)  cos* . 
Since  v,0mv,,  i4#"t^  +  (Siv#ni/2iy,cos2*  +  (Bl«;,n</Sj 
xcos#(t/}  -  (B,v,n</2flJ’  •in,#),/*,  and  80*  8  -  fit  -  # 


*8-0/,  and  thus  x0-  (v10/0)sin80*x  -  (t>x/0)sin8, 

♦  *  (ux/0)  sinflj*  (8-  0/),  we  then  obtain 

/,(v,  r,  t)=f,(vt,e,v„x,t) 

=  (m/2t7',),/,exp(-(m/27',){t^  +  vj 
+  (Btw.n//2B,)*  cos2*  +  (B,b.O f/B,) 
x  cos*  [vf  -  (B,b40//2B,)*  sin**  ]l/ *})  (14) 

Tx(/)  =  T,[l  +  (efc£,/mw)*f*/8].  (15) 

(ii)  Second  harmonic  cyclotron  resonance  heating 
(n  =  2):  We  keep  only  the  term  in  (9)  and  (10)  and  let 
Jl(w)mw/2  for  lwl<l,  (wa2vl/c).  Since  cos*  =  (q 
kp)/w,  and  sin#  =  ka/w,  we  can  then  solve  (9)  and  (10) 
and  obtain  q  +  W(/)=q[coshUa/)-  cos*  slnhGtof)]  and 
fca(/)  =  qsin*  sinh(kat),  where  a  =  Btv<(/4B(.  We  now 
have  v,0=v,,  o’0=  t>^[cosh(2fea/)  -  cos#  sinh  (2  fej /)]"*, 
a*  Blv,/4B0  and  ♦  =  *(*-  (ux/0)sin8)  +  2(8-  Of).  Thus, 
we  have 

/,(v,  r,  /)  8,  v.,x,  t)  =  (m/2»7j),/* 

x  exp(-  (m/2T,){i£  +  v*[cosh(2fea/) 

-  cos*  sinh(2fea/)]'1}) ,  (16) 

TtU)  =  T,  expl(  V2m)(fe*eE,/mu>*)*f*] .  (17) 

The  results  show  that  for  long  interaction  times,  the 
second  harmonic  heating,  which  increases  exponentially 
in  f*,  is  more  efficient  than  the  fundamental.  Note  that 
equation  (15)  is  different  from  that  given  in  Ref.  11 
where  TL  increases  as  /. 
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ABSTRACT 


Cyclotron  resonance  heating  of  bouncing  electrons  by  an  obliquely 
incident  wave  field  is  analyzed.  Continuous  interaction  between  elec¬ 
trons  and  wave  fields  throughout  the  entire  range  of  electron  excursion 
has  been  considered  in  the  analysis.  The  results  show  that  the  bouncing 
notion  of  electrons  serves  to  alleviate  the  detuning  effect  of  frequency 

mismatch  and  efficient  heating  can  be  achieved.  Kinetic  temperature  is 

2 

found  to  Increase  algebraically  in  t  for  the  fundamental  resonance  and 
exponentially  in  t  for  the  second  harmonic  heating.  Due  to  the  interac¬ 
tion  between  electrons  and  the  parallel  component  of  the  wave  electric 
field,  heated  electrons  tend  to  focus  to  the  midplane  and  their  excur¬ 
sion  amplitudes  also  become  fllamented  in  the  steady  state. 


I.  Introduction 


Interactions  between  high  frequency  electromagnetic  fields  and 
plasmas  have  been  studied  extensively  in  connection  with  wave  propaga¬ 
tion  in  the  ionosphere  and  with  the  studies  of  plasma  heating  with  high- 
power  microwaves.  The  commonly  recognized  mechanism  for  achieving  effect¬ 
ive  electron  heating  is  the  electron  cyclotron  resonance  heating  (ECRH) . 
Recently,  ECRH  has  been  adapted  to  toroidal  devices  and  received  consid¬ 
erable  attention  due  to  the  encouraging  results  reported  for  experiments 

12  3  4 

in  the  Elmo  Bumpy  Torus  (EBT)  *  and  in  Tokamaks  ’  . 

Several  theories  have  been  developed  to  explain  the  observed  heating 
rates.  Usually  the  problem  can  be  approached  by  two  different  ways.  One 
approach  is  to  study  the  wave  propagation  in  plane-stratified  plasma. 
Techniques  of  geometrical  optics  have  been  employed  for  this  purpose, 
and  it  was  found  that  full  wave  theory  is  required  near  the  resonant  sur¬ 
faces**’  Another  approach  is  by  solving  the  equations  of  motion  of  a 

8—15 

single  electron  moving  in  the  wave  fields  and  averaging  the  results 
over  an  Initial  distribution.  However,  for  studying  the  localized  heat¬ 
ing  phenomenon,  for  instance  the  heating  of  electron  rings  in  the  EBT  de¬ 
vice,  the  latter  approach  may  give  better  insight  of  the  physics.  It  is 
Important  to  note  that  the  heating  efficiency  depends  on  the  rate  of  net 
energy  gain  throughout  the  entire  period  of  particle-wave  interaction.  To  make 
this  point  clear,  we  recall  that  the  power  delivered  from  the  wave  to  the 
charged  particle  is  the  dot  product  of  the  electric  field  force  jf  the 
wave  acting  on  the  charge  particle  and  the  velocity  of  the  charge  particle. 
Even  if  the  force  and  the  velocity  are  initially  large  and  in  the  same 
direction  (so  giving  large  initial  heating  rate),  if  the  phase  angle  differ¬ 
ence  between  these  two  vector  increases  linearly  in  time,  the  dot  product 
will  oscillate  becoming  zero  on  the  average.  Thus,  the  heating  efficiency 
will  be  zero  on  the  average.  This  is  why  resonance  heating  is  effective, 
because  at  exact  resonance  the  phase  angle  between  heating  force  and  particle 
velocity  remains  constant  (fundamental  resonance  case)  in  the  linear  region 
and  continuous  energy  flow  from  wave  to  charge  particle  (or  vice  versa  can 


be  achieved.  However,  an  exact  resonance  condition  is  impractical  in 
reality,  and  a  small  mismatch  to  the  resonance  condition  means  a  sttong 
deterioration  in  the  "effective"  heating  time  and  so  in  heating  efficient, 
the  detuning  effect. 

There  are  two  interlinked  sources  of  detuning:  one  is  the  inhomogenity 
of  the  background  magnetic  field  and  the  other  is  the  motion  of  the  part¬ 
icle  guiding  center. 

In  a  uniform  magnetic  field  case,  for  Instance,  if  the  heating  wave 

is  not  exactly  normal  incidence  (i.e.  k^  i  0)  and  also  if  there  is  no 

other  imposed  force  field  in  the  system,  the  k^z  phase  angle  in  the  wave 

fields  can  introduce  a  detuning  effect  of  frequency  mismatch  to  the  charged 

particles  and  the  continuous  heating  scheme  breaks  down.  Unless  other 
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effects  such  as  collisions  and  stochasticity  *  are  set  into  the  system 
to  randomize  the  detuning  phase,  effective  heating  can  not  be  achieved. 

In  the  analyses  of  such  cases  one  normally  assumes  that  a  single  part¬ 
icle  will  repeatedly  pass  through  a  narrow  resonance  zone  and  gaining 

a  finite  but  small  velocity  Av.  After  passing  the  resonance  zone  the  part¬ 
icle  drifts  through  a  large  region  of  detuning,  suffering  a  large  phase 
shift  A0  without  gaining  energy.  After  many  such  passings  through  the 
resonance  zone,  the  scatter  of  the  A 9’s  is  presumably  large,  that  an  iso¬ 
tropic  outward  diffusion  in  the  v^  space,  and  hence  a  net  gain  of  the  rms 
value  of  velocity,  will  result.  Heating  is  thus  achieved. 

The  purpose  of  this  paper  is  to  present  a  different  view  point.  We 
believe  that  the  continuous  particle-wave  interaction  throughout  the  en¬ 
tire  range  of  the  particle  excursion  is  important.  Both  the  magnitude 
and  the  phase  of  the  velocity  evolve  continuously  whether  the  particle  is  in 
or  out  of  the  resonance  zone.  Net  energy  gain  can  be  fundamentally  more 
efficient  than  that  of  the  stochasticity  approach.  To  bring  out  the 
physics  of  this  continuous  interaction  we  choose  a  model  for  which  a  self- 
consistent  analysis  is  possible.  We  assume  a  uniform  background  field 
but  the  particles  are  confined  in  the  direction  parallel  to  the  magnetic 


field  by  an  .externally  Imposed  field  to  form  a  simple  harmonic  os¬ 
cillator.  Thus  the  detuning  effect  appears  in  the  k^z  term  of  the 
phase  function.  In  Section  II  we  first  formulate  the  problem  and  de¬ 
rived  the  coupled  equations  for  the  slow  varying  amplitude  and  phase  of 
both  the  transverse  and  parallel  motions.  In  Section  III,  the  results 
are  analyzed  in  detail,  giving  the  heating  rates  and  specific  features 
on  the  parallel  bounce  motion.  In  the  final  section,  a  discussion  on 
the  main  features  of  the  physics  involved  is  presented. 


II.  Heating  Rate  Equations 


The  equations  of  motion  of  electrons  are 


_d_ 

dt 


(1) 


^v— Jct+i  vx(B  +  zBo)]-u»gZ  z  (2) 

2 

where  -*•»  z  Is  a  harmonic  bouncing  force  acting  on  electrons  and  introduced 
®  -*■  -*■ 

by  an  external  force  field,  and  G  and  B  are  the  transverse  wave  fields  given 
by 

E“(zk^  -  x  (E/k)cos(l^x  +  k^^z  -  uit) 

4.  “ 

B  ■  -y (kcE/u) cos (k^  x  +  k^z  -  ut) 

In  order  to  provide  some  physical  insight  into  the  mathematical  pro¬ 
cedure  employed  to  study  the  problem,  we  first  decompose  Eq.  (2)  into  trans¬ 
verse  and  parallel  components  and  write  the  transverse  equation  formally 
in  integral  form: 

-i(0  -<flt)  -i<0  -fnc)£  i<0  *)t’> 

vx  -  i  vy  -  v^e  °  +  (n/uj)(kcE/Bo)e  °  ^  dt'  e 

°  (3) 

[(lyrfA-v  Jco.Oc.V  +k11z,-U)t') 


-5- 


The  integrand  of  (3)  represents  the  instantaneous  Lorentz  force 
due  to  the  heating  wave  seen  by  the  gyrating  electrons  at(x',  z'),and 
at  time  t*.  Because  of  the  gyrating  motion,  x*  is  thus  an  oscillatory 
function,  nearly  sinusoidal  at  frequency  0,  with  slowly  varying  ampli¬ 
tude  and  phase.  Therefore  the  cosine  function  in  (3)  may  be  formally 
expanded  into  an  infinite  series,  each  term  successfully  smaller  and 
oscillating  at  the  Dopier  shifted  frequency  m-mO.  If  the  heating  wave 
frequency  is  one  of  the  harmonics  of  the  cyclotron  frequency  (w  =  nfi) , 
the  m*n  term  is  the  slowest  term  and  dominates  over  all  other  terms 
which,  after  Integration,  are  smaller  by  at  least  a  factor  (fit)  ^  and 
can  be  neglected.  If  further  the  effect  of  k^z  could  be  neglected, 
e.g.  k^*0,  this  term  would  be  non-oscillatory  (refer  to  as  the  secular 
term).  Exact  resonance  is  achieved  and  heating  will  result.  However, 
if  k^z  cannot  be  ignored  such  as  the  case  of  free  streaming  particles, 
the  phase  k^v^t  gives  rise  to  a  fast  oscillating  force  resulting  in  de¬ 
tuning.  The  main  theme  of  this  paper  is  to  deal  with  the  effect  of  a 
bouncing  motion  (contrasting  to  free  streaming)  and  to  show  that,  in 
this  case,  efficient  heating  can  be  achieved.  One  may  employ  a  similar 
argument  used  in  dealing  with  the  gyrating  motion  to  deal  with  a  sinu¬ 
soidal  oscillation  in  z.  One  can  then  define  a  set  of  self-consistent 
resonance  trajectories  having  the  form 


x  "  x0  -^o^  sin  9o  +  (vjL/n)  sin(0Q  +  f  +  flt) 
vx  *=  vx  cos(eo  +  t  +  fi  t) 

vy  “  Xi  sin(6Q  +  f  +  fit) 

z  -  Z(t)  Sin($  +  u)„t)  (5 

... 

where  u  is  the  bouncing  frequency  and  the  slow  time  varying  functions  v  (t), 
D  1 

T(t),  Z(t)  and  ♦(t)  are  governed  by  the  equations,  determined  in  a  self- 
consistent  way  by  substituting  (5)  back  into  Eqs.  (3)  and  (4).  The  details 
of  this  derivation  are  shown  in  the  Appendix. 
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dt  vj.  "  ^eE/m)J0(kj1Z)[nJn(k^vi/n)3/(k^vL/f3)]cos  (nip  +  $q) 

7~t  ♦  “  -(kllkx/k)(eE/^)Jo(kii2>Cjn^va/n)/(kava/n)]sin(n^*  ) 


ftZ  -  -|cj<»’  -  Z2)%  |j1(ku2)|/Z 

#(t)  ■  +  (vio/mB)(Z^  -  Z2>*3/2^) 


with  the  initial  conditions  v^(o)  *  v^,  $(o)  *  0,  Zq  =  (z2  +  v^2/u>2)1'^2,  and 
«(o)  *  tan  1^“>bz0/vz0)  respectively,  and  C^**  (k^/k^kXcE/B^  (0/u>B)  J^k^v^ADsin  ^ 


w*iere  *  \  *o  ”  eo  +  neo'  an<*  subscripts  "o'  represents  initial 

quantities. 

III.  Analysis 

Equations  (6)  and  (7)  can  be  combined  to  obtain  an  invariant  relation 
Jn(kiVl/n5  s  in  (nip  +  4>q)  ■  I  *  constant  in  time 

Eq.  (10)  shows  that  the  result  of  resonance  is  to  change  v^  and  ^contin¬ 
uously  and  simultaneously,  following  a  closed  contour  in  the  polar  v^-n^ 
coordinate  space.  However,  the  part  of  the  contour  farthest  from  the  origin 
reaches  into  the  region  for  too  large  for  the  relativisitc  effect  to  be 
neglected.  Hence,  within  the  non-relativistic  regime,  the  heating  is  a  con¬ 
tinuous  process  as  will  be  demonstrated  and  further  discussed  later.  One 
notices  that  here  in  (10)  v^  (and  thus  Z)  does  not  appear  explicitly.  This 
Is  because  when  the  axial  excursion  becomes  bounded  due  to  bouncing,  the  time 
varying  phase  difference  "k^z"  between  wave  field  and  the  gyration  motion 
of  electrons  also  becomes  a  bounded  function.  After  the  Bessel  function  ex¬ 
pansion,  this  detuning  effect  only  appears  as  a  factor  of  JQ(kj^Z)  *n  re” 
ducing  the  strength  of  resonance  coupling  {for  example,  see  Eq.  (6)),  but  not 
giving  a  detuning  frequency  any  more.  This  modification  in  the  coupling 
strength.  Instead  of  mismatch  frequency,  will  cancel  but  when  taking  the 


(  -"-i&rtfia ,.1,'jaftr. 


t 


racio  of  Eq.  (6)  to  Eq.  (7)  for  deriving  the  adiabatic  invariant  relation. 
Therefore,  the  resultant  invariant  relation  (10)  does  not  involve  v^;  and 


v^  and  v^  only  couple  to  each  other  weakly  so  that  the  change  rate  of  v^ 

depends  on  v  through  J  (k,,Z)  factor.  Another  adiabatic  invariant  rela- 
z  o  11 

tion  for  Z  can  also  be  derived  by  Eqs.  (6)  and (8).  Since  there  is  no  need 
to  have  this  relation  to  obtain  the  final  results,  this  relation  is  not 
presented. 


Physically,  one  can  describe  the  effect  of  the  bouncing  motion  to  the 
resonance  heating  in  the  following  way.  The  power  delivered  from  the  wave 
to  the  transverse  energy  of  electrons  is  determined  by  the  dot  product  of 
and  the  transverse  Lorentz  force  F^  .  In  the  gyrating  frame,  one  may 
assume  v^  is  more  or  less  fixed  in  one  direction  but  Fj^  moves  around 
v.  due  to  the  time  varying  phase  difference.  If  the  phase  difference  is 
a  linearly  increasing  (or  decreasing)  time  function,  then  F^  is  rotating 
around  v^  at  the  detuning  frequency.  During  half  period  electron  is  gain¬ 
ing  energy  but  giving  the  absorbed  energy  back  to  the  wave  during  the  next 
half  period,  and  thus  the  net  power  delivered  from  wave  to  electrons  is  zero 
on  the  average.  The  total  energy  gained  by  electrons  is  then  limited  to 
the  fake  heating  which  is  inversely  proportional  to  the  square  of  the  detun¬ 
ing  frequency,  a  very  small  value.  Without  collisions  or  stochastic  effect 
etc.,  heating  cannot  be  continuous.  However,  if  electrons  are  bouncing  in 
the  axial  direction,  thus  k^z  becomes  a  bounded  function.  In  other  words, 
be  vibrating  about  v^  instead  of  gyrating  around  v^  .  Therefore, 
is  non-zero  over  time  averaged.  There  is  a  continuous  net  power 
flow  from  wave  to  electrons  or  vice  versa,  depending  on  the  initial  condi¬ 
tions.  But  as  shown  later  electrons  will  always  gain  energy  continuously 
when  averaged  over  the  initial  conditions.  Continuous  heating  is  achieved. 
This  is  shown  by  the  perpendicular  heating  rate  equation  derived  from  Eqs. 
(6)  and  (10)  as 


v  «U1 


d  2 
dt  vj. 


2(k11(U/k^k)(eE/m)Jo(kuZ)CJ^(k^vL/n)  -  I2]* 


(ID 
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Equation  (11)  should  be  considered  simultaneously  with  Eq.  (8)  and  can 

be  solved  numerically.  However  the  excursion  amplitude  Z(t)  evolves  very 

slowly  in  time,  Eq.  (11)  alone  may  be  used  to  determine  the  heating  rate 

corresponding  to  the  initial  phase  of  the  heating  by  letting  Z«Zq.  If  one 

further  recognizes  that  the  argument  of  the  Bessel  function  (k  v  /Q)  is 

0(v/c)  and  thus  J  (x)  s  -=j-(£) ,  Eq.  (11)  can  then  be  integrated  analytically, 
n  n!  l 

For  the  cases  of  n  *  1  (fundamental)  and  n*2  (second  harmonic)  we  obtain: 

For  n  *  1 : 


2  2 


v  -  v  +  (eE/m)(k  /k)J  (k  Z  )cos  6  v 
1  1  O  '  11  olio  oil 


For  n  ••  2 : 


+  (eE/2m)2(k11/k)V(k11Zo)t2 


2  2 

\  “\o  Icosh(bt^  +  cos  sinh(bt)  ] 


(13) 


where  b  *>  (eE/2mft)  (k,  ,k.  /k)  J  (k..Z  ). 

11  1  o  11  o 


The  above  two  equations  give  the  time  evolution  of  the  perpendicular  velocity 
v^  of  a  typical  particle  with  initial  conditions  V^Q»  (imbedded  in  4>q) 
and  an  initial  excursion  amplitude  ZQ. 

To  obtain  the  kinetic  perpendicular  temperature  T^(t)  (12)  and  (13) 
must  be  averaged  over  an  initial  distribution  of  the  particles.  In  the  per¬ 
pendicular  direction  the  averaging  is  straight  forward  since  the  initial  con¬ 
ditions  appear  only  as  a  phase  4>q.  With  the  averaging  in  the  parallel  direc¬ 
tion  formally  expressed  we  have: 


n  -  1 


Tl(t)  ‘  Tlo  +  <Jo<»t11Zo)>(e2E2/8m)  (k11/k)2t2 


(14) 


and  n  ■  2 


T( (t)  ■  T,  <cosh  bt> 
1  lo 


(15) 


In  the  parallel  direction,  the  particles  have  both  kinetic  and  potential 

energies  (harmonic  oscillators) .  The  excursion  Z  is  thus  directly  related 

2  2  ° 

to  the  total  particle  energy:  mu  Z  /2.  In  a  Boltzmann  distribution,  we 

o  O 

have 


1 

oz  J 


Jo(kllZo)exp 


Z  dZ 
o  o 


exp[-  (k2n1oz/ w2m)  ]  IQ  (k^T^m) 


k.,  2irT 
11  oz 


) 


1/2 


(16) 


where  I  (x)  is  the  modified  Bessel  function,  and  T  is  the  initial  parallel 
o  oz 

temperature. 

With  this  result,  (14)  becomes 


T.(t)-T,  +1.58  f_XT-1/2(k../k)P  t2  keV  (17) 

A  lo  B  oz  11 

where 

fg  :  bouncing  frequency  in  KHz 

X  :  wave  length  of  heating  wave  in  meters 

T  :  parallel  temperature  in  electron  volts 
zo  2 

P  :  heating  wave  power  density  in  watts/100  cm 

t  :  time  is  ps 
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For  the  second  harmonic  case  (n“2),  the  averaging  of  cosh  bt  is  not 

easy.  However  because  of  the  Jensen  inequality, it  can  be  shown  that 
2  1/2 

<cosh  bt>>cosh(<b  >  t].  Thus 

2  1/2 

T.(t)>T.  cosh[<b  >  '  t)  (18) 

A  -  A  O 

<b2>1/2-0.099(k,/k)(k  /k)1/2(f  AP)1/2(T  )'1/4  (ws)-1 

A  11  B  ZO 

The  units  of  various  quantities  are  the  same  as  those  given  in  (17) . 

The  results  show  that  for  long  interaction  time,  the  second  harmonic  heating, 
which  increases  faster  than  an  exponential  increase  in  t,  is  more  efficient 
than  the  fundamental,  a  similar  conclusion  to  the  previous  results  obtained 
under  the  situation  of  exact  resonance  but  without  a  bounce  field.  It  thus 
implies  that  the  bounce  motion  of  electrons  in  the  axial  direction  can  indeed 
alleviate  the  detuning  effect  of  frequency  mismatch  as  mentioned  before. 

"o  study  the  nonlinear  evolution  of  the  excursion  amplitude,  Eq.  (8) 

should  be  integrated.  We  first  consider  the  case  |k^z|  <  1,  and  let 

J  (k,  ,Z)  =1  and  J,  (k, ,  Z)  ;  k,  ,Z/2,  thus  Eq.  (8)  is  solved  to  give 
o  11  1  11  11 

Z2(t)  =Z2  cos2(|c  |k„t/2)  (19) 

o  O  11 

This  result  shows  that,  for  electrons  with  initially  small  excursions,  the 

excursion  amplitude  itself  will  oscillate  slowly.  For  arbitrary  values  of 

the  initial  excursion  Z  ,  Eq.  (8)  is  solved  numerically  and  the  results  are 

o 

shown  in  Fig.  1.  Basically,  for  k^Zo<3.84,  the  first  finite  zero  of  J^(x), 
the  oscillatory  behavior  of  Z(t),  Eq.  (19),  is  seen.  For  larger  values  of 
Zq,  k^Z(t)  decreases  monotonically  and  approaches  the  next  zero  of  J^(x) 


asymptotically.  Since  the  successive  zeros  of  J^(x)  are  separated  by 
about  it,  the  bounce  amplitude  of  each  particle  is  thus  decreased  by  an 
amount  less  than  X^/2,  representing  a  small  reduction  of  the  parallel 
temperature. 

The  steady  state  bounce  amplitude  Z(“)  is  seen  to  become  discretized 

to  the  zeros  of  J, (x).  For  reasonably  large  X,,  and  Z  the  discretization 
1  11  o 

of  Z(")  will  show  up  as  heavy  concentration  of  particles  at  the  turning 

points  where  the  instantaneous  velocity  is  zero.  This  seems  to  be  in 

agreement  with  experimental  observation^.  In  the  present  model  of  a  simple 

harmonic  oscillation  in  z,  it  can  be  shown  that  the  density  is  mildly  sing- 
-1/2 

ular  as  (z-z.)  .  However  a  slight  spread  of  w  will  destroy  the  sing- 

1  D 

ularity  yielding  a  smooth  density  peak.  For  X^  small  and  with  an  in- 
ltally  cool  plasma,  the  separation  between  the  discretized  values  of  Z(<*>)  is  small. 
A  concentration  of  density  at  the  mid-plane  will  then  appear  since  all 
particles  pass  through  the  mid-plane  and  more  particles  have  small  excur¬ 
sions.  Physically,  this  is  due  to  the  interaction  of  gyrating  and  bounc¬ 
ing  electrons  with  the  wave  field.  The  phase  relation  between  the  motion 
and  field  evolves  nonlinearly  so  that  the  axial  secular  force  on  the  elec¬ 
trons  point  inward  to  the  mid-plane. 
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IV.  Discussion  and  Conclusion 


Eqs.  (17)  and  (18)  represent  rather  efficient  rates  of  heating.  These 

results,  although  obtained  through  the  approximation  J  (x)  *  (x/2)n/n! ,  are 

n 

valid  before  reaching  the  relativistic  regime.  A  discussion  of  the  physics 
background  is  thus  in  order.  In  case  of  fundamental  heating,  we  can  express 
(12)  on  the  simultaneous  evolution  of  v^(t)  and  <t>(t)  by  the  triangular  dia¬ 
gram  of  Fig.  2,  where  $0  is  determined  entirely  by  the  initial  conditions 
of  the  particle.  It  is  seen  that  all  particles  will  gain  in  v^(t)  eventu¬ 
ally  even  for  the  most  unfavorable  initial  condition  cos  <f>o  =  —  1  .  -i(»(t) 
evolves  to  approach  asymptotically , in  inphase  with  the  wave  field.  For 
the  case  of  the  second  harmonic  heating,  n  ■  we  rewrite  (13)  as 


2.  .  1  2  .  .  .  bt  ,  .  -bt, 

(t)=-rv.  [(1  +  cos  <J>  )e  +  (1  -  cos  4>  )e  ] 

1  /  lo  o  o 


(20) 


It  is  seen  again  that  all  particles,  except  the  set  of  zero  measure 
cos  «J>0  "=  +1  (for  b>o  respectively),  eventually  will  gain  in  v^(t).  The 
effect  of  all  particles  gaining  v^  differs  fundamentally  from  the  stochastic 
approach  where  the  particles  diffuse  in  the  v  -  ^  plane.  For  fundamental 
heating,  temperature  increases  in  t  to  a  power  less  than  2. 


Finally,  we  would  like  to  refer  to  the  very  important  work  by  Jaeger 

14  15 

et  al  and  by  Lieberman  and  Lichtenberg  on  ECH  in  a  parabolic  mirror. 

A  comprehensive  study  was  reported  in  the  two  papers  cited  which  included, 

among  other  aspects  of  ECH,  an  extensive  amount  of  numerical  work  from  which 

most  of  quantitative  conclusions  were  drawn.  The  detuning  effect  treated 

was  that  due  to  the  inhomogeneous  magnetic  field  of  a  mirror  and  the  k^z 

effect  was  not  considered.  In  spite  of  the  difference,  their  numerical 

2  2 

solution  also  indicates  a  t  growth  of  <v^> .  However,  it  is  difficult  to 
deduce  from  the  paper  the  rate  of  heating  for  comparison.  In  contrast  to 
the  continuous  interaction  approach  employed  here;  the  impulse  approximation 
was  used.  Whether  continuous  interaction  can  indeed  predict  similar  im¬ 
proved  heating  rate  in  a  mirror  still  remains  an  open  question. 


Since  the  background  field  configuration  for  practical  devices  is  usually 

■ore  complicated  than  the  model  used  in  the  present  analysis,  only  qualita¬ 
tive  explanation  to  the  experimental  observation  is  expected.  For  a  quanti¬ 
tative  comparison,  a  more  elaborate  model  should  be  adapted. 
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•  FIGURE  CAPTION: 


Figure  1 


Figure  2: 


:  Tine  Evolution  of  the  Excursion  Amplitude  for  cases 
different  Initial  values,  where 

|Co|  -  (k^1/kLk)(cE/Bo)(G/aiB)Jn(kLvio/fi)sin 


Evolution  of  v^-phasor  during  Fundamental  heating. 


Appendix 


To  derive  the  equations  for  the  slow  time  varying  functions  v^(t) , 
♦(t),  Z(t)  and  *(t),  we  substitute  the  resonance  trajectory  (5)  into  Eqs. 
(3)  and  (4).  The  Integrand  of  the  integral  on  the  right  hand  side  of  Eq. 
(3)  can  be  expressed  by  Bessel  function  expansion  as 


[(kj^u/k  )  -  vz]cos(k^x  +  k11z -ut) 
i  l  <klf*2-l'  (kn2)  Jt  l(k  v  70) 

I0>  I 


i(£-n) (6q  +  fit)  i[(£-l)iR$  1  i£’(W-oj  t) 
e  e  °  e  B 


+  Ji+l(klvl/n)e 


-i(£-n)  (0  +I2t)  -i[(i+l)<H-4>  ]  -il^frh^t)! 


where  4>Q  *=  k jX^kjV^o/fiisin  $o  +  n0g»  vz  *  WgZcos($  +  u^t)  ,u  *  nlJ,  and  the 

expansion  e  3  -  ][  J  (a)e*P^  is  used.  Therefore,  the  slow  time 

P  «  -  «  P 

varying  part  of  (Al)  can  be  obtained  by  setting  l-  n  and  £'  =o,  the  result  is 


i(6  +Jlt) 

O 


[(k1iM/k2)-vi]cos(kA*  +  kii«-  wt)| 


slow  part 


-  (kn“/k2)J0(knZ)j^~-  Jn(kJvi/n)cos(n*  +  $o> 
+  i  J^(k^v^ /n)sin(ni|i  +  $D)je  ^ 


Substituting  (A2)  Into  (3) ,  yields 


v  e'^-v  +  (k. .  /k)  (eE/ra)  /*dt '  J  (k  -Z') 

1  |0  11  O  O  11 

cos(n$'+$  )  +  i  J'  (k.  v' /Jl)sin(n4»'  +  ♦  ) 
o  nil  o 


(nfl/k^v^ > Jn  ^ kivi /R) 

-i\|»' 

e 


(A3) 


from  which  after  first  time  derivative  Equations  (6)  and  (.7)  are  obtained. 

Similarly,  the  first  term  on  the  right  hand  side  of  Eq.  (4)  can  be  expanded 

as 

[v^- (k^u/k^lcosCk^x  +  k^  z-wt) 

4  I  [(W/kj^)  «(liu)/k2)]JJlCkLvi/n)Jt,(k11Z) 

1,1' 

i  i(i-n)(e  +nt)  i<£*  +  *  )  ii'(»  +  u  t)  » 

je  °  °  e  B  +  c.c  (A4) 


Thus,  the  resonance  term  is  obtained  by  keeping  only  terms  in  (A4)  with 
£  *  n  and  £*  ■  +1.  The  result  is 

j[vx-(kiU/k2)]cos(kAX+ku  z  - ut)|  resonance 

(A5) 

-  -  2  (k2  co/k  k2)J  (k  v  /ll)J.  (k.  -Z)sin(n<|i  +  4  )sin($  +  u>Rt) 
n  1  nij.  a  ii  °  D 

With  the  aid  of  (A5) ,  Eq.  (4)  becomes 
A2  7 

-~z  +  u„z»-2ujbC  J1(k..Z)sin(*  +  w_t)  (A6) 

j.  a  D  B  O  1  XI  o 


2 

where  C  ■  (k*  /k  k) (cE/B  )(n/u_)J  (k  v  /n)sin*  and  the  invariant  rela- 
o  11  1  o  Bnllo  To 

tion  derived  from  Eqs.  (6)  and  (7)  has  been  used  to  reduce  the  expression 

of  C^.  Eq.  (A6)  describes  a  harmonic  oscillation  driven  by  a  resonant 

force,  the  solution  of  it  is  given  by 

z-zqcos  u>Bt+ (vzo/u)B)sin  Ugt  +  A(t)cos (*  +  Wgt)  (A7) 

£  a  *'2 
where  A(t)  =  C^^dt'  Jo(k^Z’),  and  I#!  «  oig  and  |A/A|  «  oig  are  assumed. 

From  the  original  description  z  -  Z(t)sin(4>  +  w  t)  and  (A7),  two  self-con- 

B 

slstent  equations  for  Z(t)  and  4(t)  are  obtained 

z  +  A  cos  ♦=  Z  sin  ♦ 
o 

(A8) 

(v  /a  )  -  A  sin  ♦«=  Z  cos  4 

ZO  fi 

Uhich  are  solved  to  obtain 

L 

Z2-z^+  (v20/ug)  +  A2  +  2A[zq  cos  *-  (vzo/wB)sin  ♦  ]  =  Z2  -  A2  (A9) 

and 

♦  -  sin"1  {[zqZ  +  vz{A/«i>b ] /Z2  -  s in"1  ( [ zqZ  +  (vz0/w^(Z2  -  Z2 ) 1 / 2 ] /Z2 }  (A10) 

where  the  relation  zq  cos  ♦-  (v20/uB)sin^  “  -  A  which  can  be  derived  by  (A8) 
is  used,  Zq- (*2  +  v2o/«J)1/2,  and  A  - +(Z2  -  Z2)1/2  is  given  by  (A9) . 

Taking  first  time  derivative  on  (A9) ,  yields 


.(All) 


The  usage  of  plus  or  minus  sign  in  (A10)  and  (All)  depends  on  the  initial 
conditions.  However,  |z|  is  always  decreasing  from  Zq  [i.e.  |z|  < Zq  as 
shown  by  (A9)],  it  is  thus  convenient  to  express  (All)  and  (A10)  as  Eqs. 
(8)  and  (9),  respectively. 
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Consideration  of  the  effect  of  a  uniform  pump  field  on  the  particle  orbits  in  a  Vlasov  plasma  leads  to  a 
modified  diffusion  coefficient.  When  the  particles  oscillate  in  the  pump  field,  the  turbulent  wave  phase 
velocity  seen  by  the  particles  is  Doppler  shifted  by  multiples  of  td^fk  (a>,  is  the  pump  frequency).  Hence, 
strong  interactions  between  the  particles  and  various  components  of  the  turbulent  field  will  take  place.  It  it 
shown  that  when  the  pump  field  is  sufficiently  strong:  £,  £  (4<rn,7', )"’,  bulk  heating  can  dominate  over  tail 
heating  and  excitation  of  electrons  to  energy  levels  higher  than  the  normal  suprathermal  values  is  possible. 
This  field  strength  is  within  the  range  of  laser  fusion. 


I.  INTRODUCTION 

During  the  past  several  years,  substantial  interest 
has  been  shown  in  the  problem  of  anomalous  heating1'4 
in  a  plasma  driven  by  a  large  field  oscillating  near 
to,,.  Parametric  instabilities  are  excited  and  believed 
to  be  responsible  for  the  observed  heating  of  electron 
plasmas  in  laser  fusion,’  ionospheric  modification,* 
computer  simulation, ’■’•*  and  laboratory  experiments. 
Heating  of  suprathermal  tail  electrons  at  about  the  phase 
velocity  of  the  electron  plasma  wave  is  predicted  in 
various  theories.’ •*  Although  the  heating  of  suprather¬ 
mal  electrons  is  indeed  observed  in  experiments,7'* 
bulk  heating* '7'*  and  the  production*  of  electrons  with 
velocities  several  times  the  phase  velocity  of  the  elec¬ 
tron  plasma  wave  have  also  been  observed  in  the  labor¬ 
atory  experiments  and  computer  simulations.  Kaw 
el  al .*  analyzed  the  coupling  between  fast  plasma  waves 
and  slow  plasma  waves  through  an  independently  excited 
ion  density  fluctuation.  Bulk  heating  was  also  predic¬ 
ted.  In  this  work,  we  have  properly  taken  into  account 
the  pump  field  effect  on  the  particle  orbits,  which  has 
been  neglected  in  previous  theories.  All  wave  compo¬ 
nents  present  are  either  parametrically  excited  in¬ 
stabilities  or  driven  through  wave-wave  interactions. 
They  possess  wavelengths  of  about  the  same  order.  A 
modified  diffusion  coefficient  is  derived  in  Sec.  II.  This 
diffusion  coefficient  is  used  to  study  both  bulk  heating 
and  tail  heating.  Substantial  bulk  heating  and  heating 
at  higher  than  the  suprathermal  levels  are  predicted. 
The  ratio  of  the  energy  deposition  rates  to  the  bulk 
electrons  and  the  suprathermal  tail  electrons  is  estima¬ 
ted  in  Sec.  III.  A  discussion  is  given  in  Sec.  IV. 

II.  FORMULATION 


E(r0>  t)  =  E„,,  exp{i[k-  r„  -  (w  +m*>0)t] } 

is  the  excited  turbulent  field  with  w=ui(k),  the  frequency 
of  the  excited  ion  wave.  In  most  of  the  parametric 
theories,10  only  those  waves  with  n  =0,±1  which  cor¬ 
respond  to  the  ion  wave  and  the  up  and  down-shifted 
electron  plasma  waves,  respectively,  are  considered. 
However,  as  the  pump  level  is  increased,  higher  fre¬ 
quency  electron  fluctuations  (m  +  l)u>0±  w,  m  -  1,2,. ..  , 
can  also  be  excited  through  the  coupling  of  the  pump 
field  and  those  fundamental  instabilities  (n  -0,il  modes). 
Since  these  higher  frequency  fluctuations  are,  in  gen¬ 
eral,  off-resonance  modes,  the  thresholds  for  the 
parametric  excitation  of  these  higher -frequency  fluc¬ 
tuations  are  high.  Therefore,  these  fluctuations  are 
not  instabilities  but  driven  waves  sustained  by  the 
pump  field  and  the  fundamental  instabilities.  The  rela¬ 
tions  between  these  higher  frequency  fluctuations  and 
these  fundamental  instabilities  are  derived  in  the 
Appendix.  Following  the  procedures  by  Dawson  and 
Oberman,11  and  by  Kuo  and  Cheo,1’  we  now  substitute 
the  relations 


T«=T  sint*V.  T'*r -%2$cosw0t 

m  O<*>o  m  aw0 

into  Eq.  (1),  and  define 

e(r,  /)  =E[r  -  (qaK0/mow20)  cosio0f ,  f] 


€i.»exp{i[k.r -(cu  +  fu>0)/]},  (2b) 

i  * 


where 


We  consider  the  problem  of  turbulent  heating  of 
plasmas  by  a  strong  electromagnetic  pump  wave  through 
parametric  excitation.  The  Vlasov  equation  can  be  ex¬ 
pressed  as 


(I) 


where  a  ^electron,  ion;  E,(f )  =  E, cosu>0f  is  a  prescribed 
spatially  uniform  pump  field,  and 


Thus,  Eq.  (1)  becomes 

(i7+T,fe+^e(r-/),^)/(Tir’,)'0'  (4> 

where  /(v,  r,  I )  */„[ t  +(<?oE0/m<,u>0)  sinu>0f ,  r  -  (qJEe/ 
m0b>2)cos<a>0f,(],  the  distribution  function  in  the  (r  ,») 
system. 

We  now  proceed  to  derive  the  renormalization  quasi- 
linear  diffusion  equation.  We  notice  that  the  pump  field 
E,  is  removed  from  the  transformed  Vlasov  equation  (4). 
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We  assume  that  each  component  of  the  turbulent  field  all  other  effects  which  are  to  be  ignored.  Expanding  the 

<  carries  a  random  initial  phase  such  that  the  ensemble  spatial  dependence  of  /(1)  in  a  Fourier  series 

average  over  the  initial  phases  =0  for  all  k  and  _ 

integers  f.  The  distribution  function  is  written  as  the  =  £  /( 4exp{t[k.r  -  (w  ♦  fw0)f)} , 

sum  of  four  parts  as  done  by  Dupree11  1  •* 

/»{/)+/<*>  ♦/„  =  (/)♦/“’  ♦/„,  (5)  and  substituting  (2),  (5),  and  (6)  into  (4),  wet 

where/11’ =/**’♦/,(/)  is  the  average  distribution  tunc-  /a  b\...  Jls.  8 

tion,  /'*’  is  the  phase  coherent  response  to  the  electric  \at+y’ar)'''~~ma  by'  * 

field,  /  describes  the  phase  incoherent  portion  due  to 

the  nonlinear  wave-wave  interaction,  and  /.describes  with  /, .  */{  ‘1  +/,  » ,  where 


fu\r,r,t)  =  £  /,  >  exp{’[it«*-  -  (w  ♦  iw0)/]}  , 

I  .* 

and  substituting  (2),  (5),  and  (6)  into  (4),  we  obtain 


with  /M  */{'j  *fi  *  •  where 


rn„  i[k.  v  -  (<j  +  Ud0)\  By  'J ’  m\  fri,  «!*•  v- (w  + 1  w0)j 

Ed  t  _ _ l  d  .(c) 

By  '  **’  "*'  «[(k-k')-  v  -(w-  w'  +  (l-n)w0)l  ’  By'1* 

r>  _ 1 _ ^  _9_  _ _ f  _L /(e)* 

+  ma  i[k.v-(ui +  !<*{,)]  2-«av  ’  *’*'*'•  ij(k-k').  v  -  (to  -  w  +  (l-n)u>0)l  ’  By'  ’ 


»?.  «(i*v  -(w+fu>0)J  ?  8r  + 


Equation  (8)  can  be  solved  iteratively,  and  f\‘j,  can 
be  expressed  in  terms  of  </)  including  all  the  higher 
order  terms.  The  result  of  the  expansion  can  be  ob¬ 
tained  by  solving  an  integro-differential  equation  for 
the  averaged  propagators  in  a  way  similar  to  that  of 
Dyson.14  We  express  approximately 

yhJa.Jx  f  rfr  exp{-i[k.  v-(w +  /a>0)]r -5- k-0  •kr’} 

*  9ft  _  J 


*«,,•£</>•  (10] 

If  only  wave-particle  resonant  interactions  are  consid¬ 
ered,  we  may  substitute  (10)  into  (7)  and  neglect  the 
contributions  from  the  nonlinear  wave-wave  interac¬ 
tion  terms  contained  in  (9).  Thus,  the  renormalized 
quasi-linear  diffusion  equation  is  obtained  as 


X  dr  £ 


(£*t-8-7)</>'£-0.-£</>. 


xexpf-tfk.v-fw+foJoJlr-dJk-  D0.kr*}  (12) 

is  the  renormalized  diffusion  tensor  which  includes 
resonance  broadening  due  to  the  turbulent  field  «  and 
the  effect  due  to  the  presence  of  the  coherent  pump  field 

From  Eq.  (3),  we  write 

mtn' 

x  ^1  -»'(a«K^»a  ®»'a)  •  (13) 

As  shown  in  the  appendix,  waves  at  frequencies 
(m«-l)u0tiii,  where  m  =1,2,. .. ,  are  sustained  by  the 
coupling  of  the  pump  field  and  parametrically  excited 
instabilities  at  frequencies  u>  and  w„±u>.  Therefore, 
the  ensemble  average  of  (13)  can  be  expressed  in  terms 
of  the  Intensities  of  those  instabilities  (n=0,±l  com¬ 
ponents)  as 


Substituting  (14)  into  (12),  the  diffusion  coefficient  of 
the  turbulent  spectrum  of  the  parametric  instabilities 
is  derived.  This  diffusion  coefficient  is  different  from 
the  one  derived  in  Refs.  IS  and  16,  in  that  the  summa¬ 
tions  of  the  Bessel  functions  appear  here.  In  the  pre¬ 
viously  obtained  diffusion  coefficient,  only  tail  heating 
by  I  ♦»!.»!*  is  predicted.  However,  when  this  effect  of 


the  coherent  pump  field  is  included,  additional  coupling 
between  waves  and  particles  is  also  obtained.  This 
can  be  observed  from  expression  (14).  For  instance, 
the  f>0  term  corresponds  to  the  part  of  the  diffusion 
coefficient  for  the  bulk  heating  of  electron  plasma. 

This  term  includes  contributions  from  both  ion  waves 
and  electron  waves.  The  |J|  *  1  terms  contribute  to 
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heating  o{  suprathermal  electrons  as  shown  in  the  work 
of  Katz  et  al.*  Also  predicted  here  is  the  production  of 
electrons  of  energies  at  higher  than  those  of  the  supra- 
thermal  electron  levels  given  by  the  |f  |  =  1  terms.  The 
production  of  those  higher  level  electrons  is  also  ob¬ 
served  in  the  computer  simulation  work  of  Kruer  and 
Dawson.'  The  additional  coupling  can  be  understood 
from  the  fact  that  particles  are  oscillating  in  the  pump 
field;  wave-particle  interactions  are  thus  Doppler 
shifted  by  the  multiples  of  the  oscillating  frequency. 
This  is  similar  to  the  higher  harmonics  in  cyclotron 


resonances.  If  the  resonance  broadening  effect  is  neg¬ 
lected  in  (12),  the  modified  quasi-linear  diffusion  co¬ 
efficient  is  then 

O-^E  <e.>£i*..>6[k-»-(u)  ♦/«,)].  (11 

mo  I  .» 

Again,  we  observe  additional  spikes  of  the  diffusion 
coefficient  in  velocity  space,  not  given  in  previous 
theories.  If  only  up  to  the  first  order  of  the  Bessel 
functions  are  kept,  and  if  only  the  diffusion  of  elec¬ 
trons  is  of  interest  (|a,|«  |a,|),  (14)  becomes 


<«!.*« +Jj(a,)|0J^|’) , 

{<  -i.»<  -i .*)a  ♦-1^1*1 » 

«.  ..«/.*>  -  I **  [i  -  Js  A(o.)  (i  -  ■/«<.>) I  ♦,.*!* . 

<«-a a> - kk[l  - ^’.(<0 (l  -  •'!<«.>)  V.) I  <#»-i j*  I*- 


In  reality,  the  exact  diffusion  coefficient  given  by  (12) 
would  include  peaks  at  different  points  in  velocity  space 
each  broadened  to  a  finite  interval  around  the  peak. 
Similar  to  (15)  these  peaks  are  widely  separated  at  a 
distance  |u>0/*l  which  has  a  value  about  several  elec¬ 
tron  thermal  speeds.  Therefore,  these  intervals  hardly 
overlap  each  other,  and  we  may  then  decompose  the 
diffusion  coefficient  D  into 

0  =  0, +0, , 


l  dT 

x  exp[-i(k-T- w)t  -(i)k-  0,-kr'J, 

°,=E  °... 


/  drE  E  <«i 


x  exp{-t[k •  r  -  (u>  +fu)0)]r  - (j)k.  D„.kr'}.  (18) 

0,  and  D,  are  the  diffusion  coefficients  for  the  supra¬ 
thermal  tail  and  for  the  bulk  heating,  respectively.  The 
significance  of  O,  as  shown  in  (17)  is  that  |^,ia>|',  the 
intensity  of  the  electron  plasma  waves,  which  usually 
carries  more  energy,  do  contribute  to  bulk  heating.  No 
previous  theories  have  demonstrated  this  effect.  For 
oscillating  two-stream  instabilities,  they  are  always 
the  dominant  terms  for  bulk  heating.  In  decay  instabil¬ 
ities,  they  become  dominant  when  </'(a,)>  (uj/w0)J\(o,) 
as  given  by  the  Manley-Howe  relation.  Moreover, 

I also  contributes  to  the  production  of  higher 
energy  electrons.  This  effect  is  included  in  Eqs.  (I6d) 
and  (16e).  Physically,  the  instability  fields  can 
again  couple  to  the  pump  field  in  order  to  sustain  fields 
oscillating  at  2wt±u  and  thus  transfer  energy  to  elec¬ 
trons  at  the  higher  energy  level.  However,  more  im¬ 
portant  is  that  the  time  variation  of  the  fields  >,lt*  ex¬ 
perienced  by  the  oscillating  electrons  in  the  pump  field 


is  not  at  a  single  frequency  as  would  be  seen  by  a  sta¬ 
tionary  observer.  Instead,  it  contains  a  Doppler  shifted 
spectrum.  Thus,  the  resonant  interaction  between 
electrons  with  speeds  at  about  (2u,tu)/li  range  and 
p,,,,  also  occur.  Although  contributions  from  these  two 
mechanisms  are  out-of-phase,  they  never  cancel  each 
other.  The  Doppler  shift  effect  always  dominates.  The 
arguments  given  are  also  true  for  the  |J|>2  terms. 

Thus,  significant  heating  and  production  of  energetic 
electrons  is  also  predicted  as  the  pump  power  is  in¬ 
creased.  The  heating  rate  and  hence  the  ratio  of  the 
energy  deposition  rates  to  the  bulk  electrons  and  the 
suprathermal  tail  electrons  will  be  analyzed  in  the 
next  section. 

III.  ELECTRON  HEATING  RATE 

In  order  to  simplify  the  analysis,  we  only  consider  the 
one-dimensional  case.  In  a  homogeneous  plasma,  the 
spectra  of  parametric  instabilities  are  symmetric  in 
the  dipole  approximation,  and  hence  the  diffusion  coef¬ 
ficient.  Thus,  the  average  background  (/)  in  the  os¬ 
cillating  frame  is  also  symmetric  in  velocity  space  if 
the  initial  distribution  is  assumed  to  be  a  Maxwellian 
without  drift.  The  zeroth  moment  of  (II)  gives  the  con¬ 
tinuity  equation,  and  the  first  moment  of  (11)  indicates 
that  there  will  be  no  drift  velocity  generated  by  this 
symmetric  turbulent  spectrum.  The  second  moment  of 
(11)  thus  determines  the  rate  of  energy  gain  due  to  the 
turbulent  heating  as 

/  <J)j^lvD(v)]dv ,  (19) 

where  c  «(m/2)  j u*(/)di;.  We  now  follow  Dupree’s 
argument,"  the  integrals  of  (17)  and  (18)  may  be  ap¬ 
proximated  by  step  functions  Cf(x): 
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where 

/W3»'»  .  VM  r.  .  /»e» 3‘”  .  \»» 

D°°  Urn*#*  4“°/  ’  ° *  “  \4m*^‘  *"■/  ’ 

u-EkMi*.  <—=53 u«..-r + ic— >ia). 

*  * 

*0  is  the  central  wave  number  of  the  spectrum.  Upon 
simplification  of  the  diffusion  coefficient,  Eq.  (19)  be- 


F1G.  1.  Initial  background  distribution. 


srw« 


■"  I  -“'■))  '{’{-i.  *“■*))]  -fe  *"'•)[(/ 

-  -■}>]-(=?  -■#(=?  -"'■))]} 


8  8 

'I7£»+I7£r. 


and 

^-£,-m  /  </>|>^)rfo  +  2m£D.[(^  -W ^  -  (=g>  ♦**)  (/(^^r,))]  •  (24) 

r  “ 

By  viewing  Fig.  1,  the  right-hand  sides  of  (23)  and  (24)  Maxwellian.  Equations  (25)  and  (26)  may  be  expressed 

are  proportional  to  the  areas  of  the  shaded  regions  I  explicitly  to  be 

and  II,  respectively.  Hence,  (23)  and  (24)  may  be  ap-  .  r  /  /  \\ 

proximated  a, 

— -  2mD.(^  [('  (*;  -  “'•))  -{/  (7,  *  "'•))]  .  vd,{(A0))  (/in'.)))!  u(w, 

-.{(/<«» -l[(.  (£-»,))  ,2m0„U.fc)-|L//fc)) 

•('Ml'M).  -  -&***+■-& 

and 

and 

(28)  —  A/«b\\  (28) 

To  calculate  the  ratio  of  the  energy  deposition  rates  ‘V*#/  ”  T-*  Y  \*o//  ’  ' 

into  the  bulk  electrons  and  the  suprathermal  tail  elec-  where  ^  . {m/2tTyn  exp(_mt WtJ  is  assumed,  and 
irons,  one  has  to  derive  (/)  from  Eq.  (1 1).  Although  „  >  j  terro,  (26)  are  neg,ected.  Therefore, 

the  diffusion  coefficient  is  assumed  to  be  constant  in 

each  heating  region,  as  time  evolves,  the  distribution  /i9c,M.\/9cTH,.\~l^/Ji(g.)yjru»c  / 

function  (/)  will  no  longer  be  Maxwellian.  The  govern-  \  8 1  /[  9 1  )  w0(a,)/  (L 

ing  equation  (1 1)  may  be  solved  numerically.  However,  idW'.UlW.  -  w/k.)"]  *)  /  mw*  \ 

we  present  a  reasonable  estimate  of  the  distribution  of  ♦  1  1 — I  *****  ( cxp(2jfcly  )  • 

the  heating  energy,  by  deriving  the  initial  energy  depo-  '  0  *° 

sition  rates  in  the  following.  We  assume  that  (/),.,  is  (29) 
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At  moderate  pump  power  levels,  i.e., 
n0  =  E0/(4»nor.0),/,<l1  I  a.  |«  1 . 
and  Ji(a#)a  aj'l  and  J0(a,)=  1.  Thus, 


. 

is  the  total  energy  in  all  the  plasma  waves  and  Pmy/vtl 
is  the  {actor  of  order  one  or  larger  by  which  the  pump 
intensity  exceeds  its  threshold  value.”  Substituting  (30) 
into  (29)  we  obtain 

(sr)(iir)"-(^){[“I*,",(55r 

*  (!«*)  V.  -  •  <’■> 

As  an  example  we  show  that  bulk  heating  can  become 
equal  to  tail  heating  even  at  moderate  pump  power 
levels.  We  consider  the  case  of  oscillating  two-stream 
instabilities:  gj  =0  and  u>„=  w„.  Using  we 

obtain 


&!“)  (£««.)"  *l- 


when  ij0* 0.187  P‘in  * 0.187.  Using  ho/fe4,=0.2,  we  ob¬ 
tain 


when  qo=0.62  P*1/,<0.62.  However  at 
*•/*'•  •0M- 

when  tj0»  1.52  P“tn  «  1.52,  this  value  is  in  the  strong 
pump  power  region.  Since  strong  Landau  damping 
occurs  in  the  region  ka/ku  *  0.25,  instabilities  have 
maximum  growth  rate  in  the  region  k/kt,*0.2.  There¬ 
fore,  from  the  numerical  values  given  here,  it  is  rea¬ 
sonable  to  expect  bulk  heating  to  dominate  over  tail 
heating  for  a  sufficiently  strong  pump:  |£9|»(4im<,TM),/*. 

IV.  DISCUSSION 


frequency  and  k  is  the  wavenumber  ol  the  turbulent 
field.  Therefore,  the  parametrically  excited  electron 
plasma  waves  do  interact  with  those  bulk  electrons 
effectively.  Moreover,  the  production  of  electrons  at 
several  phase  speed  of  the  electron  plasma  waves  is 
also  predicted.  Numerical  estimates  show  that  bulk 
heating  can  dominate  over  tail  heating  for  a  sufficiently 
strong  pump. 

In  the  present  analysis,  it  is  seen  that  the  value  of 
a,  =  (k/fe<4)[ £0/(4 rnaT.)ln  1  is  of  critical  importance  to 
which  part  of  the  velocity  spectrum  be  heated.  At 
small  values  of  «„  bulk  heating  can  be  larger  than  tail 
heating  provided  that  £0/(4*n0T,)l,J  a  1.  Bulk  heating 
is  maximum  at  a,  =  1.84.  At  higher  values  of  a,  the 
higher  order  Bessel  functions  can  no  longer  be  neglected 
and  heating  to  much  higher  velocities  than  those  of  the 
suprathermal  electrons  are  possible.  For  £„( 4en0T,)',/’ 
=  1  the  intensity  of  the  pump  is  about  10u  W/ctn1  for 
glass  laser  and  10”  W/cm’  for  CO,  laser  (7,  ■  10  keV). 
This  is  within  the  presently  available  intensities  of  both 
glass  and  CO,  lasers. 
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APPENDIX 

We  now  derive  the  relations  between  those  higher 
frequency  fluctuations  [(m  +  l)u>0±u>,m  =  l,2,...]  and 
the  fundamental  parametric  instabilities  (h>,<ua±u))  by 
using  Vlasov’s  equation  and  Poisson’s  equation.  The 
relations  are  only  derived  at  the  linear  level,  namely, 
we  linearize  Eq.  (4);  thus,  the  density  perturbation  6fa 
in  linear  response  to  the  total  turbulent  field  c  in  the 
oscillating  frame  is  governed  by 

(ji •  i)  6/»(v’ r* ■  it <(r  ■ f) •  &«{v)  -  lA1) 

where  f^v)  is  an  unperturbed  background  distribution 
in  the  oscillating  frame,  it  may  be  simply  assumed  to 
be  Maxwellian.  If  we  define  an  unperturbed  trajectory: 

( d/dt)r  =v  and  (d/df)v=  0,  and  thus  v=v'  and  r-r' 
=▼(1-0,  where  r=r(f),  r' =r(f’)  and  v  =v(f),  v'  =  v(f'); 
hence,  (Al)  can  be  integrated  and  the  result  is 

Til  A 

x  exp{fk»  r  -  i[u>  +(n  +m)w0]/} 


The  perturbed  orbit  plasma-turbulence  theories  of 
Dupree”  and  of  Weinstock"  predict  the  broadening 
of  the  diffusion  region  as  the  turbulence  level  increases. 
Hence,  the  wave-particle  interactions  become  more 
effective  than  the  result  of  conventional  quasl-llnear 
diffusion.  If  the  turbulence  is  parametrically  excited 
by  a  coherent  pump,  the  particle  orbit  is  perturbed  by 
the  pump  field  in  addition  to  the  turbulent  fields.  The 
oscillation  of  the  particle  due  to  the  pump  field  will  also 
broaden  the  interaction  between  the  turbulent  fields. 

This  is  due  to  the  fact  that  the  phase  velocity  of  the  tur¬ 
bulent  waves  as  viewed  by  the  particles  is  Doppler  shif¬ 
ted  by  integral  multiples  of  <•>,/*,  where  ui0  is  the  pump 


.  E.Jt.8/8v/MI 

i(k  •  v  -  (o>  +(n  +  m  )<*>„])  ‘ 


Since  only  electrostatic  fluctuations  are  considered,  we 
may  define  E,.,  “-ik^.^  and  substitute  it  and  (A2)  into 
Poisson’s  equation,  which  yields 

[»•?* 


[»•?* 

-  -  £ 4°  £  <«' --*-.<«.)  Jm.,  («„) 

a  *  l 

m 

. /w+0»+m-f)wo\  . 

*W{ 
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where 


With  the  aid  of  (A4)  and  (13),  Eq.  (14)  is  obtained. 


f"  7ZTZW 


and 


*1,  =4 M„«Vrs  . 


Since  those  high-frequency  waves  are  strongly  off- 
resonance  to  the  plasma  modes,  we  may  simply  as¬ 
sume  that  they  are  sustained  by  the  fundamental  in¬ 
stabilities  Instead  of  being  excited  parametrically. 
Hence,  they  have  a  much  smaller  amplitude  than  the 
fundamental  Instabilities  and  (A3)  may  be  approximated 
as 


and 


0i».* 


x 


»-i 


J-.(a,)d..t-i(o.) 

(2n  —  p  -  l)1  u$ 


♦i.» 


for  n  *2,3 


(A4) 
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